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Abstract 

We study the amount of interference in random quantum algorithms using a recently derived 
quantitative measure of interference. To this end we introduce two random circuit ensembles com- 
posed of random sequences of quantum gates from a universal set, mimicking quantum algorithms 
in the quantum circuit representation. We show numerically that these ensembles converge to the 
well-known circular unitary ensemble (CUE) for general complex quantum algorithms, and to the 
Haar orthogonal ensemble (HOE) for real quantum algorithms. We provide exact analytical formu- 
las for the average and typical interference in the circular ensembles, and show that for sufficiently 
large numbers of qubits a random quantum algorithm uses with probability close to one an amount 
of interference approximately equal to the dimension of the Hilbert space. As a by-product, we 
offer a new way of efficiently constructing random operators from the Haar measures of CUE or 
HOE in a high dimensional Hilbert space using universal sets of quantum gates. 
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I. INTRODUCTION 

It is generally acknowledged [ij that quantum information processing differs from classical 
information processing fundamentally in two ways: the use of quantum entanglement and 
the use of interference. While quantum entanglement is undeniably 2J of crucial importance 
in tasks like quantum teleportation [3j , and has evolved into a scientific field of its own (see 
4| for a recent review), its role in quantum algorithms is less clear. Large amounts of 
entanglement are necessarily produced in any quantum algorithm that provides a speed-up 
over its classical analogue, but it remains to be seen if the entanglement is a by-product 
rather than the fundamental basis of the quantum speed-up [5J. 

Interference on the other hand has received comparatively little attention in the context 
of quantum information processing. It has been used for a long time to test the coherence 
of quantum mechanical propagation 0, EJs], and it has been proposed as a tool to cre- 
ate entanglement between distant atoms 9(, but its role in complexity theory is virtually 
unexplored [101 ] . 

In contrast to entanglement, which is a property of quantum states, interference charac- 
terizes the prop agation of states. A quantitative measure of interference was introduced very 
recently in ll|. It was shown that a Hadamard gate creates one basic (logarithmic) unit of 
interference (an "i— bit" ) . Basically all known useful quantum algorithms, including Shor's 



and Grover's algorithms 12|, [l3j start off with massive interference by applying Hadamard 



gates to all qubits. However, the two algorithms differ substantially in the amount of in- 
terference used in their remaining non-generic part: while the factoring algorithm uses an 
exponential amount of interference also for that part (in fact a number of i-bits close to the 
number of qubits), only about 3 i-bits suffice for the rest of the search algorithm, and that 
number is asymptotically independent of the number of qubits. 

The existence of an interference measure makes it meaningful for the first time to ask 
the following questions: How much interference is there typically in a quantum algorithm 
running on n qubits? How is the interference distributed in an ensemble of quantum algo- 
rithms? What is the average interference, what its variance? Are these values different if 
the algorithm has a real representation? 



2 



II. INTERFERENCE DISTRIBUTION IN THE CIRCULAR RANDOM MATRIX 
ENSEMBLES 



In order to talk about the statistics of interference, the ensemble needs to be specified. 
It is well known that any quantum algorithm (i.e. any given unitary transformation in 
the tensor product Hilbert space C 2 ®") can be approximated with arbitrary precision by 



a „ of q _ gates act,,, on at .nost two qub ,ts at the t im e Q M . More 
precisely, a universal set of quantum gates is formed by a fixed U(A) transformation, such as 
the controlled-NOT gate (CNOT) acting on two arbitrary qubits, in conjunction with the 
set of all U(2) transformations of any single qubit. Alternatively, any quantum algorithm 
may be represented by only real (i.e. orthogonal) matrices at the price of doubling the size 
of the Hilbert space, with a universal set of quantum gates consisting of the Hadamard 
gate and the Toffoli gate [l?], Q- Without any further prior knowledge of the quantum 
algorithm, it is natural to chose algorithms from Dyson's circular unitary ensemble (CUE) for 
unitary algorithms, and from the so-called Haar orthogonal ensemble (HOE) for algorithms 
representable by an orthogonal matrix [3]. CUE corresponds to an ensemble of unitary 
matrices which is flat with respect to the Haar measure dn^iU) of the unitary group U(N); 
and HOE to an ensemble of orthogonal matrices which is flat with respect to the Haar 
measure cI[in(0) of the orthogonal group O(N), where g?/^v(0) is invariant under right and 
left orthogonal transformations (dfi(0) = d^{ViOV2) for any two orthogonal matrices V\ 
and V2). We will provide numerical evidence further below that CUE and HOE represent 
more realistic quantum circuits indeed very well, once the number of quantum gates is large 
enough. 

The measure of interference introduced in 11[ reduces in the case of unitary propagation 



by a iV x iV matrix U with matrix elements Un~ in the computational basis to 

N 

1(U)=N-J2 \Uik\\ (1) 

i,k=l 

with < T(U) < N — 1. Of the two characteristics of interference, coherence and superposi- 
tion of a large number of basis states (equipartition) , only the latter distinguishes different 
entirely coherent quantum algorithms representable by a unitary matrix. The maximum 
amount of interference is reached for any quantum algorithm which spreads out each com- 
putational basis state equally over all computational basis states, whereas the interference 



is zero for a mere permutation of the computational basis states. 

We have numerically calculated the distribution of interference Pcue,n(Z) of N x N 



matrices from CU 
unitary matrices 



u sing the Hurvitz parametrization for creating large ensembles of random 
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20] . Figure [T] shows the result for N between 2 and 8. With growing 
N, the distribution becomes increasingly peaked on a value close to N. For N = 2 the 
distribution can be easily calculated analytically. We parametrize U2 with four angles a, ip, 
X chosen randomly and uniformly from the [0, 2%[ and cp = arcsin(£ 1//2 ) with £ random and 
uniform from [0,1], 

(cos (pe 1 ^ sin pe lx \ 
•v ■ J • (2) 

— sin pe 1X cos ipe 1V I 

Thus, I(U 2 ) = 4(£-£ 2 ), and 

PcnE , 2iI) = -j±=, (3) 

in very good agreement with the numerical result. Figfj] indicates that for sufficiently large 
N all quantum algorithms will typically contain the same amount of interference of order 
N. This is confirmed by an exact analytical calculation of the two lowest moments of 

n 

the interference distribution. Invariant integration over the unitary group [21[ gives closed 
formulas for integrals of the type 

2' m i 1 7T \^ m ^\TT \^ m 3 



Z U!N ( mi ,m 2 ,m 3 ) = J dti N (U)\U ilh \^\U ilh r*\U h32] 

mi!m 2 !m 3 !(^-2)!(iV- 1)! 
(N + m 1 -2)\(N + m 3 -2)! 
(N + mi + m 3 - 2)! 
(N + mi + rri2 + m 3 — 1)! ' 

where d/j,N(U) is normalized to f d^^iU) = 1, and 11,12,31,32 are arbitrary indices. This 
leads to the average interference 

(T) U>N = N - N 2 Z U>N (2, 0, 0) = N(l - ^- T ) (4) 
N -^?N-2, (5) 
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with (. . .)u,n — f dfJ,N{U){- ■ •)• The 2nd moment can be found from 
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) U;N = {N{N-l)) 2 Z UiN {2,0,2) 

+ 2iV 2 (iV- 1)Z(2,2,0) 

+ N 2 Z UtN (A,0,0) 
N 2 + 2N-1 
(JV+l)(JV + 3)" ( } 

Thus, the standard deviation of the interference distribution in the CUE ensemble, 



vanishes like ~ 2/iV for large A/". 

Figure [1] also shows the interference distribution Phoe.n{Z) for the HOE ensemble, rel- 
evant for quantum algorithms representable with purely real (orthogonal) N x N matrices. 
We constructed this ensemble numerically by diagonalizing real symmetric matrices drawn 



from the Gaussian orthogonal ensemble (GOE) 22], which for the relatively small matrix 

nn 

sizes turned out to be more efficient than Hurvitz's method [19J, 120] . Remarkable is the 
symmetric structure of the interference distribution for N = 2, whose analytical form is 
easily obtained from 2x2 rotation matrices with uniformly distributed rotation angles, 

For N > 4, the distribution becomes mono-nodal, and more and more peaked with increasing 



N . The method of invariant integration can be generalized to the HOE ensemble 23|. The 
result corresponding to eq.Q reads 



Z , N (mi,m 2 ,m 3 ) = dfi N (0)(O lin ) mi (O ilj2 ) m2 (O i2j2 ) 



in?, 



2 2 - A T(^)r(Hf^)r(Hf*i) 
r(jv- i)r( JV+mi + m3 ~ 1 ) 

* N+mi+m2+m :i ^ 

where m 1 ,m 2 ,m 3 are all even, T means Euler's gamma function, and dfijsi{0) is normalized 
to J dfi^(0) = 1. The average interference in the HOE ensemble is then given by 

(1)o,n = N - N 2 ZoA4, 0, 0) = N(l - 

N -^N-3, (9) 




FIG. 1: (Color online) Distribution of interference in the CUE ensemble (left) and HOE ensemble 
(right) for N = 2 to N = 8 in steps of 1 (curves from left to right). The full lines for N = 2 
represent the analytical results, eqs.(|3]) and (jHJ). 

with (. . .)o,n — f dfi>N(0)(- ■ •)• Thus, a real quantum algorithm of the same size N, drawn 
from HOE, contains on the average asymptotically slightly less interference than a unitary 
one drawn from CUE. However, since the size of the Hilbert space has to be doubled to 
express an arbitrary complex algorithm as a real one [18], about twice as much interference 
is needed to run the real algorithm. The second moment 

2 

(iV(iV-l)) 2 Z OiiV (4,0,4) 



N 



leads to the variance 



0,N 



+2iV 2 (iV-l)Z o ^(4,4,0) 

+N 2 Z o ,n(8,0,0) 
3N(-A + 3N(N + 5)) 

(N + 1)(N + 2)(N + 6) ' 
2AN(N - 1) 



(10) 



(11) 



(N + 2) 2 (N 2 + 7N + 6) 
of the interference for the HOE ensemble. Thus, the standard deviation cto,v(2^) decays as 
~ 2y / 6/iV for large N, and therefore practically all algorithms drawn from HOE contain for 
large N an amount of interference X ~ iV — 3. 
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III. RANDOM CIRCUIT ENSEMBLES 



Shor's algorithm was recently shown to lead to CUE level statistics [241 ] . whereas the 
quantum Fourier transform alone has a regular spectrum (it is a fourth root of the identity 
matrix), and so does Grover's algorithm, which is to good approximation a sixth root of the 



identity matrix [251 ] . It is therefore natural to ask to what extent are the CUE and the HOE 
ensembles representative of realistic quantum algorithms? 

To answer the above question, we introduce two random quantum circuit ensembles, 
the random unitary circuit ensemble (UCE), and the random orthogonal circuit ensemble 
(OCE), constructed to resemble realistic quantum algorithms with n g randomly chosen gates 
as follows: 

• For quantum gate number i (1 < i < n g ), decide whether to apply a one-qubit gate 
(with probability p) or a multi-qubit gate (with probability 1 — p). 

• If gate i is a one-qubit gate, chose randomly, uniformly over all qubits, and indepen- 
dently from all other gates the qubit on which the gate is to act, and pick as gate a 
random unitary 2x2 matrix from CUE for the construction of an UCE algorithm, or 
the Hadamard gate for building an OCE algorithm. 

• If gate i is a multi-qubit gate, chose randomly, uniformly over all qubits, and indepen- 
dently from all other gates a control qubit (two control qubits) and a target qubit, and 
apply the CNOT gate (the Toffoli gate) to these qubits for UCE (OCE), respectively. 

• Repeat this procedure for all gates i = 1, . . . ,n g and concatenate the obtained gates 
to form the entire quantum algorithm. 



A similar ensemble of random quantum circuits was introduced in [26(, where, however, 
one random constant depth gate was iterated, and the entangling gate was constructed 
from simultaneous nearest neighbor interactions. Nevertheless, according to [27| , one might 
expect at least an exponential convergence to CUE (HOE) also for UCE (OCE), respectively, 
and this is what we are going to show numerically. 

We first examine the convergence by comparing the distribution P(s) of nearest neighbor 
spacings s of the eigenphases <pi of the N x N unitary matrices. For large N, and average 
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s normalized to unity, CUE leads to a P(s) well approximated by the Wigner surmise 22] 



32s 
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Pw(s) = ^ r e- is2 ^. (12) 



Deviations are of order 10 



281 ] . For n = 4, the minimum number of gates that leads to 
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an approximately constant density of eigenphases, such that unfolding the spectrum 
is unnecessary, is n g ~ 10. For even smaller numbers of gates strong peaks at <p = and 
(p — TV arise in the density of states corresponding to a predominance of real eigenvalues, but 
otherwise the density is already flat. Fig. [2] shows P(s) for UCE for n = 4 and several values 
of n g (n r = 10 5 realizations). For small n g , P(s) has a strong peak at s — 0. The rest of 
the distribution is between the Poisson result of uncorrelated phases, P(s) = exp(— s), and 
the Wigner surmise Pyy(s). The peak at s = becomes smaller and smaller as the number 
of gates increases, and at the same time a more and more pronounced maximum at s — 1 
arises, resulting in a distribution which rapidly approaches the Wigner surmise, eq. (fT2l) . 
For n g = 40, P(s) is virtually indistinguishable from Pw(s). We examine the convergence 
quantitatively with the help of the quantity 

F * = J o (VPuce(s) - y/P w (s)) 2 ds = 2 (l - jf y/P(s)P w (s) ds^j (13) 

which measures a squared distance between the (square roots of) the level spacing distribu- 
tions Puce(s) of UCE and Pw(s). Fig. [3] shows F s as function of n g for UCE for various 
values of p and n = 4 qubits, obtained from 10 3 random algorithms. For p different from 
and 1, F s decays to a good approximation exponentially as ~ exp(— b(n,p)n g ), with a 
rate b that depends on p and n, before saturating at a small level largely independent of p. 
The latter is due to the numerical fluctuations in P(s) present for any finite N, as is easily 
checked by varying N. The finite precision of Pw{s) for N > 2 sets another lower bound 
on the values of F that can be possible achieved. Fig. [3] also shows that b(n,p), as obtained 
from a fit of \nF s to a linear function of n g between F 2 = 2 and F 2 = 0.1 has a maximum 
around p = 0.5. The convergence rates decrease with increasing n, and the maximum of the 
convergence rate as function of p shifts to somewhat smaller values of p. 

Numerical evidence presented in [19( indicates the same form of P(s) for HOE as for CUE, 
eq. (fT2"|) . in particular a quadratic level repulsion P(s) oc s 2 for s 1. We have examined the 
convergence of OCE to HOE based on P(s) Pw(s) as well, and have found similar results 
as in Figj3l However, it is clearly not possible to determine the limiting ensemble based 
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FIG. 2: (Color online) P (s) for n g = 10 (red), 15 (orange) 20 (yellow) and 40 (green — increasing 
values of P{X) in this order) for UCE with n = 4 and p = 0.5, n r = 10 5 matrices, compared to the 
Wigner surmise Pw( s ) (dashed black line) and the Poisson result, P(s) = exp(— s) (dashed-dotted 
black line). 




FIG. 3: (Color online) Convergence of P(s) for UCE to P\y( s ) as function of n g for n = 4 
qubits, n r = 10 3 , and different values of p: p = 0.1 black squares, p = 0.2 red squares, p = 0.3 
green diamonds, p = 0.4 blue triangles up, p = 0.5 indigo triangles left, p = 0.6 brown triangles 
down, p = 0.7 grey triangles right, p = 0.8 violet pluses, and p = 0.9 magenta Xs (left). Rate 
of convergence b as function of the probability p for n = 3 (circles), n = 4 (squares), n = 5 
(diamonds), and n = 6 (triangles up). 

on -P(s) alone. We therefore also examined directly the interference distributions P{T) for 
both random circuit ensembles, as -P(X) is, after all, what we are interested in. FigJH shows 
how the interference distribution of UCE for n = 4 evolves between n g = 10 to n g = 100 
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FIG. 4: (Color online) Interference distribution Puce,n(Z) (left), for n g = 
10,20,30,40,50,60,70,80,90,100, maxima increasing in this order, and Poce,n (right) for 
n g = 20, 50, 70, 100 (red, orange, blue, maroon, respectively) compared to the Pil) of the circular 
random matrix ensembles (black lines). All curves have n = 4 (N = 16); n r = 10 5 random 
realizations were used for UCE and OCE, n r = 10 7 for CUE and HOE. 

from a broad flat distribution to the strongly peaked interference distribution of CUE. The 
interference distribution Poce,n{Z) for OCE fluctuates much more for a given number of 
gates compared to the one for UCE, Puce,n(Z), but rapidly approaches Phoe,n(I) as well. 
To examine the convergence quantitatively, we define the quantity Fj as in f[T5j) . but with 
P(s) and Pw(s) replaced by the interference distributions Puce,n(Z) and Pcue,n(Z) for 
UCE and CUE (by PoceA 1 ) and ProeA 1 ) for 0CE and HOE). Note that P C ue(X) and 
-Phoe(^) have now to be computed numerically as well. We did so for the same dimension 
of Hilbert space N = 2 n considered for UCE and OCE. We used n r = 10 7 realizations for 
n G {4, 5}, for both CUE and HOE, as well as for n = 6, HOE; n r = 5 • 10 6 for n = 7, and 
n r = 10 6 for n = 8 (HOE); and n r = 4 • 10 5 for n = 6 (CUE). The number of realizations 
chosen for the OCE and UCE ensembles was 10 5 , with the exception of n r = 4 • 10 5 for 
n = 6, UCE. Figj5] shows the results for Fi(n g ) for n = 5. The curves for the other values 
of n examined (n G {4, 6, 7, 8}) look very similar, but the convergence slows down with 
increasing n. 

Fi{n g ) for both UCE and OCE is very well fitted by a Gaussian, at least up to the point 
where the crossover to the saturated behavior occurs. This is in contrast to the exponential 
convergence of P(s). A fit of In Fj(n g ) to a — c(n,p) n 2 g in the range 2 > Fi{n g ) > 0.01 yields 
c ~ 10 -4 for 4 < n < 6, with a maximum of c around p ~ 0.5 for both OCE and UCE (see 
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FIG. 5: (Color online) Convergence of P(I) for UCE (left) and OCE (right) to the interference 
distribution of CUE and HOE, respectively, as function of n g for n = 5 qubits (N = 32) and 
different values of p. Same symbols as in Figj3l 
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FIG. 6: (Color online) Dependence of the fitted Gaussian convergence rate c on the probability p 
for a single qubit gate (a random U(2) for UCE (left) and Hadamard for OCE (right)) for various 
numbers of qubits (circles, squares, diamonds for n = 4, 5, 6), respectively, and in addition triangles 
up and triangles down for n = 7, 8 for UCE. 

Figi]). 

The exponential (or even Gaussian) convergence of the random circuit ensembles to the 
corresponding circular ensembles as function of the number of quantum gates provides a 
new way of economically creating random unitary operators with a flat distribution with 
respect to the appropriate Haar measure in an exponentially large Hilbert space. The 
method will work on any quantum computer on which the relevant universal set of gates is 
available. The OCE is particularly interesting, as the only randomness resides in the indices 
of the qubits selected as entry of the Hadamard or Toffoli gates. About 2n g such random 
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numbers % in the range < % < n — 1 are needed for p = 1/2. This offers the possibility 
to construct truly random operators, with only a logarithmic overhead of qubits: A small 
register of n a ~ 3 log 2 n + 1 auxiliary qubits can be brought repeatedly into superposition 
of all computational states by applying a Hadamard gate to each qubit; then the register 
is measured in the computational basis and gives a random number. The outcome of one 
particular qubit (the highest significant, say), can be used to choose between the Hadamard 
and Toffoli gates, and the remaining bits specify the qubit (s) on which to act. Obviously, 
one might as well use the actual work qubits to generate these random numbers initially 
and store them for later use in the quantum algorithm. 

The question remains open whether the method presented here is efficient in the sense 
that the number of quantum gates n g needed for a given fidelity Fj or F s increases at most 
polynomially with n. This would require that the exponents c(n,p) and b(n,p) decay no 



faster than an inverse power of n for a given p. The same problem was encountered in [26|, [27 ] 
and so far no definite answer has been found. In order to address the question numerically, 
much larger values of n have to be considered. The study of the interference distributions 
is clearly not well suited for this purpose, as each algorithm (i.e. a very large dimensional 
matrix) gives only one number. 

IV. SUMMARY 

As a summary, we have introduced two ensembles of random quantum algorithms, UCE 
for general unitary algorithms, and OCE for real orthogonal algorithms. We have provided 
numerical evidence that these ensembles converge for sufficiently large numbers of gates 
and for a finite probability for both one-qubit and two-qubit gates (or three-qubit gates), 
to the well-known random matrix ensembles CUE and HOE, respectively, at least in the 
sense of coinciding level spacing distributions and interference distributions. One might 
consider these ensembles therefore as a new way of efficiently creating random unitaries 



from the corresponding Haar measure 26|, [27|. The method is universal in the sense that it 



runs on any quantum computer with a universal set of quantum gates. We have calculated 
numerical distributions of interference over the CUE and HOE ensembles, and have provided 
exact analytical formulas for the lowest moments. For large Hilbert space dimensions N, the 
interference distributions over CUE and HOE are peaked on their average values (F)u,n — 
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N — 2 and (T)o.n — N — 3, respectively, with a width that decays oc 1/N in both cases. 
Thus, randomly picked unitary quantum algorithms contain with high probability basically 
the same exponentially large amount of interference X ~ N. This result is reminiscent of 



similar findings for the amount of entanglement in a random quantum state [29( . Grover's 
search algorithm is therefore remarkably exceptional in the sense that its non-generic part 
(i.e. the part after bringing the computer into a superposition of all computational states) 
uses only a small amount of interference (the whole algorithm including the initial Hadamard 
gates produces exponential interference 
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